Second-harmonic generation in nonlinear left-handed metamaterials 
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We study the second-harmonic generation in left-handed metamaterials with a quadratic non- 
linear response. We demonstrate a novel type of the exact phase matching between the backward 
propagating wave of the fundamental frequency and the forward propagating wave of the second 
harmonics. We show that this novel parametric process can convert a surface of the left-handed 
metamaterial into an effective mirror totally reflecting the second harmonics generated by an incident 
wave. We derive and analyze theoretically the coupled-mode equations for a semi-infinite nonlinear 
metamaterial. We also study numerically the second-harmonic generation by a metamaterial slab 
of a finite thickness, and reveal the existence of multistable nonlinear effects. 



I. INTRODUCTION 



the cubic nonlinearity or its different generalizations. 
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Recent years demonstrated many advances in the de- 
sign and engineering of artificial composite structures 
with unique electromagnetic response, which have broad- 
ened significantly the range of possible wave phenomena 
not usually observed in Nature but being accessed in lab- 
oratory experiment. In particular, it has been shown 
that the composite metallic periodic structures may al- 
low realizing hypothetical and surprising materials, first 
predicted theoretically and termed left-handed media , 
which possess simultaneously negative dielectric permit- 
tivity and magnetic permeability. Indeed, the composite 
metallic structures consisting of arrays of wires and split- 
ring-resonators have been demonstrated to possess left- 
handed properties in the microwave frequency range 0. 
These composite materials received a great deal of inter- 
est due to their highly nontrivial and quite often counter- 
intuitive electromagnetic properties. In particular, such 
materials support backward waves, and they exhibit neg- 
ative refraction at interfaces 0. 

The study of linear wave propagation and linear prop- 
erties of left-handed materials is a major subject of re- 
search in this field. It is indeed the case when both mag- 
netic permeability and dielectric permittivity of the ma- 
terial do not depend on the intensity of the electromag- 
netic field. However, the future efforts creating tunable 
structures where the field intensity changes the transmis- 
sion properties of the composite structure would require 
the knowledge of nonlinear properties of such metama- 
terials, which may be quite unusual. Indeed, inserting 
small conducting elements into a periodic metallic struc- 
ture of the metamaterial brings an unique opportunity 
to enhance nonlinear response in such left-handed media 
because the microscopic electromagnetic field inside the 
structure resonant elements can become much stronger 
than an average macroscopic field. Indeed, it has already 
been shown analytically and numerically that the com- 
posite material containing nonlinear elements in the slits 
of split-ring resonators possesses a hysteresis-type nonlin- 
ear magnetic response 0, . This is an example of self- 
action of electromagnetic waves usually associated with 



Inclusion the elements with non-symmetric current- 
voltage characteristics such as diodes into the split-ring 
resonators will result in a quadratic nonlinear response of 
the metamaterial 0. This quadratic nonlinearity is re- 
sponsible for the recently analyzed parametric processes 
such as the second-harmonic generation (SHG) and 
three- wave mixing [8|. In particular, the first analysis of 
SHG from a semi-infinite left-handed medium has been 
briefly presented by Agranovich et al. , who employed 
the nonlinear optics approach and also made a state- 
ment that for metamaterials "... SHG in transmission 
is badly phase-mismatched" parametric process, and "it 
is then the SHG in reflection that is more interesting". 
(see Ref. 0, p. 165112-5). 



In this paper we consider the problem of SHG during 
the scattering from a semi-infinite left-handed medium 
(or a slab of the left-handed material of a finite extent) 
and demonstrate that the original paper by Agranovich 
et al. missed an important additional phase-matching 
condition, quite specific for the harmonic generation by 
the backward waves. With this condition, we demon- 
strate that exact phase matching between a backward 
propagating wave of the fundamental frequency (FF) and 
the forward propagating wave at the second harmonics 
(SH) is indeed possible. This novel phase matched pro- 
cess allows for creating an effective "quadratic mirror" 
that reflects the SH component generated by an incident 
FF wave. 



The paper is organized as follows. In Sec. 2 we describe 
our model including both the electric and magnetic re- 
sponses. Quadratic nonlinearity and the SHG process in 
metamaterials are discussed in Sec. 3. In Sec. 4, we 
develop the corresponding coupled-mode theory for SHG 
with backward waves and present the analysis of both 
lossy and lossless cases of this model. Section 5 is devoted 
to the case of SHG process a slab of finite-extension. At 
last, Sec. 6 concludes the paper. 
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II. MODEL 

We consider a three-dimensional composite structure 
consisting of a cubic lattice of conducting wires and split- 
ring resonators (SRR) , shown schematically in the insert 
of Fig. ^ We assume that the unit-cell size of the struc- 
ture d is much smaller then the wavelength of the propa- 
gating electromagnetic field and, for simplicity, we choose 
a single-ring geometry of the lattice of SRRs. The results 
obtained for this case are qualitatively similar to those 
obtained in more involved cases of double SRRs. This 
type of microstructured medium is known to possess the 
basic properties of left-handed metamaterials exhibiting 
negative refraction in the microwave region. 

In the effective-medium approximation, a response of 
this composite metallic structure can be described by av- 
eraged equations allowing to introduce the effective di- 
electric permittivity and effective magnetic permeability 
of the form 

eH = l-4» (1) 



fi{co) = 1 + 



Flo 



(2) 



where tu p is the effective plasma frequency, lo is a reso- 
nant frequency of the array of SRRs, F is the form-factor 
of the lattice, and to is the angular frequency of the elec- 
tromagnetic waves. The product of permittivity e and 
permeability /x defines the square of the effective refrac- 
tive index, n 2 = e.[i, and its sign determines if waves can 
(n 2 > 0) or cannot (n 2 < 0) propagate in the medium. 
Due to the medium dispersion defined by the dependen- 
cies and (J2J), the wave propagation becomes possible 
only in certain frequency domains while the waves de- 
cay for other frequencies. Metamaterial possesses left- 
handed properties when both e and [i become simultane- 
ously negative, and such a frequency domain exists in the 
model described by Eqs. Q and J3| provided co p > lo . 
In this case, the metamaterial is left-handed within the 
frequency range 



uo a < lo < min {uo p , lom} , 



lo m 



LOq 



VT^f 



(3) 



where lo p is the plasma frequency introduced in Eq. 

We assume that lom < u p , and in this case we have 
two frequency ranges where the material is transparent, 
the range where the material is left-handed (LHM), and 
the right-handed (RHM) domain for lo > lo p , where both 
permittivity and permeability are positive (shaded do- 
mains in Fig. For the frequencies outside these two 
domains, the composite material is opaque. 



III. QUADRATIC NONLINEARITY AND 
BASIC EQUATIONS 

The composite material becomes nonlinear and it pos- 
sesses a quadratic nonlinear response when, for example, 
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FIG. 1: Frequency-dependent magnetic permeability n 
(solid) and electric permittivity e (dashed) of the compos- 
ite. Two types of the regions (LHM or RHM) where the 
material is transparent are shaded. For other frequencies it 
is opaque. Characteristic frequencies loq, lom, and lo v are de- 



fined in Eqs. 
metamaterial. 



Q to 0. Inset shows the unit cell of the 



additional diodes are inserted into the SRRs of the struc- 
ture |(J , as shown schematically in the insert of Fig. ^ 
Quadratic nonlinearity is known to be responsible for var- 
ious parametric processes in nonlinear media, including 
the frequency doubling and generation of the second- 
harmonic field. In dispersive materials, and especially 
in the metamaterials with the frequency domains with 
different wave properties, the SHG process can be rather 
nontrivial because the wave at the fundamental frequency 
and the second harmonics can fall into different domains 
of the material properties. 

The most unusual harmonic generation and other para- 
metric processes are expected when one of the waves (ei- 
ther FF or SH wave) has the frequency for which the 
metamaterial becomes left-handed. The specific interest 
to this kind of parametric processes is due to the fact that 
the waves in the left-handed media are backward, i.e. the 
energy propagates in the direction opposite to that of the 
wave vector. Both phase-matching condition and nonlin- 
ear interaction of the forward and backward waves may 
become quite nontrivial, as is known from the physics 
of surface waves in plasmas @. In this paper, we are 
interested in this type of parametric wave interactions. 

In nonlinear quadratic composite metamaterials, in- 
teraction of the forward and backward waves of different 
harmonics takes place when the material is left-handed 
either for the frequency lo or the double frequency 2lo. 
Under this condition, there exist two types of the most 
interesting SHG parametric processes in metamaterials. 

Case I. The frequency of the FF wave is in the range 
loq/ 2 < lo < lom 12 and, therefore, the SH wave is gen- 
erated with the double frequency in the LHM domain 
(see Fig. QJ. For such parameters, the electromagnetic 
waves at the FF frequency are non-propagating, since 
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FIG. 2: Dispersion of plane waves k(ui) in the metamate- 
rial. Arrows show the parameters of the FF and SH waves 
corresponding to the exact spatio-temporal phase matching. 



e(w)/i(cj) < 0. As a result, the field with the frequency 
lu from this range incident on a semi-infinite left-handed 
medium will decay exponentially from the surface inside 
the metamaterial. Taking into account Eqs. and 
the depth <5 of this skin-layer can be found as 



-1/2 



< 



x_ 

17' 



(4) 



where fc|| is the tangential component of the wavevector 
of the incident wave, and A is a free space wavelength. 
For the SH wave generated in this layer, the metama- 
terial becomes transparent. In this case, a thin slab of 
a metamaterial may operate as a nonlinear left-handed 
lens that will provide an image of the source at the sec- 
ond harmonics, as suggested recently in Ref. 01 ■ This 
case has also been mentioned in the earlier paper on the 
second-harmonic generation 0. 

Case II. The FF wave is left-handed, whereas the SH 
wave is right-handed. Such a process is possible when 
ojp < 2wo (see Fig. 0. What is truly remarkable here 
is the possibility of exact phase-matching of the SHG 
parametric process, in addition to the cases discussed 
earlier in Ref. ■ The phase-matching conditions for this 
parametric process are depicted in the dispersion diagram 
of Fig. [21 for the propagating waves in the metamaterial 
where the dispersion of the plane waves is denned by the 
relation 



D(oj,k) 



0. 



(5) 



The exact phase matching takes place when 2k{uo) — 
k(2u>). Different signs of the slopes of the curves at the 
frequencies u> and 2oj indicate that one of the waves is 
forward, while the other wave is backward. 

To study the SHG process in metamaterials we con- 
sider a composite structure created by arrays of wires 



and SRRs. To generate a nonlinear quadratic response 
of the metamaterial, we assume that each SRR contains 
a diode, as depicted schematically in the inset of Fig. ^ 
The diode is described by the current-voltage depen- 
dence, 



U 

Rd 



1 



U 



(G) 



where U c and Rd are the parameters of the diode, and 
U is the voltage on the diode. Equation 10 is valid pro- 
vided U <§; U C1 and it represents two terms of the Tay- 
lor expansion series of the realistic (and more complex) 
current- voltage characteristics of the diode. 

Following the standard procedure, we consider two 
components of the electromagnetic field at the fundamen- 
tal frequency lu and its second harmonics 2ui, assuming 
that all other components are not phase matched and 
therefore they give no substantial contribution into the 
nonlinear parametric interaction. Subsequently, we write 
the general coupled-mode equations describing the simul- 
taneous propagation of two harmonics in the dispersive 
metamaterial as follows, 



AHi + e(w)/i(u)— Hi = -tnHJHa, 

cr 

uj 2 

AH 2 + 4e(2w)M(2u;WH2 = -0- 2 H 2 , 



(7) 



where the indices "1","2" denote the FF and SH fields, 
respectively; A is a Laplacian, and other parameters are 
defined as follows 



(Ti = k/2R(uj), o~2 = k/R*(uj), 
6ir (na 2 )' 



d 3 c 5 



U c R d R{uj)R{2uo) 



(8) 



where R(uj) = uJqUj 2 + iju, the asterisk stands for the 
complex conjugation, a, d are the radius of the SRRs 
and the period of the metamaterial, respectively, and 7 
is the damping coefficient of the SRR. For simplicity, we 
assume that both FF and SH waves are of the same polar- 
ization, and therefore they can be described by only one 
component of the magnetic field. In this case, Eqs. (0 
become scalar. In the derivation of Eqs. Q we take into 
account the Lorentz-Lorenz relation between the micro- 
scopic and macroscopic magnetic fields pj|. Also, it is 
assumed that the diode resistance Rd is much larger than 
the impedance of the SRR slit, i.e. Rd ^ 1/ujC 1 so that 
the resonant properties of the composite are preserved. 



IV. SEMI-INFINITE METAMATERIAL 

First, we consider a semi-infinite left-handed medium 
and the SHG process for the wave scattering at the sur- 
face. We assume that a TM-polarized FF wave is incident 
on a LH material from the vacuum, as shown schemat- 
ically in Fig. |3| Inside the metamaterial, the wave at 
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FIG. 3: Geometry of the SHG problem. Thick arrows show 
the direction of the energy flow, thin arrows - direction of 
wavevectors. Indices "i", "r", and "tr" stand for incident, 
reflected and transmitted waves, respectively. 



the fundamental frequency satisfies the dispersion rela- 
tion (JjjJ which defines the wavenumber k. As is discussed 
above, the FF wave in the left-handed medium is back- 
ward, meaning that the normal component of the wave 
vector is directed towards the surface, i.e. in the direction 
opposite to the Poynting vector. 

When the phase-matching conditions are satisfied, the 
generated SH wave has the wavevector parallel to that of 
the FF wave (see Fig. 01. However, the SH wave is for- 
ward propagating, so that the energy at this frequency 
should propagate towards the interface. When losses are 
negligible, the FF wave will be transformed completely 
into the SH wave with the energy flows in the direction 
opposite to that of the FF wave. This kind of the SHG 
process in a semi-infinite left-handed medium is charac- 
terized by two major features: (i) the efficiency of the 
SHG process may become very high, and (ii) the SH wave 
propagates in the direction opposite to that of the incom- 
ing FF wave. 



A. Coupled-mode equations 

To describe the SHG process analytically, we employ 
the coupled-mode theory and the slowly- varying envelope 
approximation for the FF and SH fields, and present the 
magnetic fields in the material in the form: 



H h2 (t,z) = a h2 (t,z)e ikl - 2Z + c.c 



(9) 



where the amplitudes of the FF and SH fields ai >2 (t,z) 
are assumed to vary slow in both space and time, i.e. 
dai t2 /dt <§; uiai t 2, and dai t2 /dz <C kai <2 . Substituting 
Eqs. © into Eqs. Q and neglecting the second-order 
derivatives, we obtain the coupled equations, 



da i 
~~dt 



dai 

dz 



ioia 1 a 2 — fldi, 



da 2 
~dt 



— - , da 2 . 2 
ST + %2^— = io 2 a x + v 2 a 2 - illa 2 , 
dz 



(10) 



where v g i <2 are the group velocities and 1^1,2 = v g i t2 Im(k) 
are linear damping coefficients of the FF and SH fields, 
respectively, 



fl = q 2 D(2u),2k)/2a 2 , 
is the phase mismatch, and 



qi = <?i 



dD(u,k) 



q 2 = o"2 



dD(2uj,2k) 



du> 



(11) 



•(12) 



The coupled-mode equations (|10|) can be presented in the 
equivalent rescaled form, 



Vgi^r- = -q\bib 2 - vibx, 
dz 

db 2 db 2 , 2 

+ v g2 -— = q 2 b x + v 2 b 2 - iilb 2 , 
oz 



~dt 



dt 



(13) 



where a x — abi, a 2 = /3b 2 , a — exp(z</<), /3 = 
exp[— i(ir/2 — 2<p)], and cj> is an arbitrary phase. 

The incoming FF backward travelling wave has the 
group velocity in the z-direction, and the phase velocity 
- in the opposite direction. The generated SH forward 
wave has both the phase and group velocities in the — z- 
direction. The FF wave propagates inside the material, 
and it loses the energy due to SHG and also due to losses 
in the medium. As a result, the SH amplitude decreases 
in the z-direction, and the boundary conditions should 
be taken in the form 61,2(00) = 0. 



B. Lossless process 

First, we neglect the metamaterial losses assuming, for 
simplicity, that v\ %2 = and Im(qi l2 ) = 0. At the large 
enough value of the phase mismatch we can neglect the 
derivatives in the second equation of the system (|13fl , and 
obtain a local coupling between the FF and SH ampli- 
tudes, 



b 2 = -l-b v 



(14) 



Substituting Eq. (|14|) into the equation for the FF field, 
we obtain a single nonlinear equation for the amplitude 
b\ with an effective cubic nonlinearity, 



dbi 



dbi 



dt + Vgl dz ' n 



101 bi 



(15) 



The general solution of Eq. 115|) can be found in the form 
bi = J 4(t,z)e l$(t ' z) , (16) 

where A(t,z) = f(t-z/v g i), <P(t,z) = (t-z/a)qiq 2 f(t- 
z/v g i)/Q(l — Vgi/a), where f(t,z) is an arbitrary func- 
tion, a is constant, and a ^ v g i. The obtained solution 
describes a stationary nonlinear wave. 
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FIG. 4: Geometry of the SHG process for a finite-width slab 
of a nonlinear quadratic metamaterial (shaded) . 



In the case of the exact phase matching, i.e. when 
£1 = 0, Eqs. 113(1 are simplified, 



dbi 



db 1 



— +V 9 1 ~d^ = - qih * h2 ' 

db 2 db 2 , 2 

at dz 



(17) 



Looking for stationary solutions of Eqs. (|17|l . we take 
d/dt = and find the integral of motion of this system, 



'12 



(18) 



Since both the waves should vanish inside the metama- 
terial, we take C = 0. Then, the solution of Eqs. (|17|) is 
found as 



9192 



1/2 



1 



1 



(2 + zo) ' 



91 (z + Zq) 



(19) 



where zo is constant. Solution (|19|l has a singularity, and 
it corresponds to an explosive instability of the paramet- 
rically interacting waves. However, this singularity does 
not appear in the corresponding boundary problem. 

Equations (fP7|l describe also the solution for the cou- 
pled FF and SH fields, where the FF wave has a localized 
profile, while the SH wave has a tanh profile, 



bt(t, z) = Aisech [{t - z/0) /T] , 
b 2 (t, z) = A 2 tanh [(< - z/0) /T] , 



(20) 




0.5 1 1.5 2 

Normalized incident field amplitude 

FIG. 5: Reflection coefficients of the generated SH wave 
(solid) and transmitted FF (dotted), as well as the transmis- 
sion coefficient of the FF wave (dashed) vs. the normalized 
field intensity. Normalized slab thickness L — 10. 



C. Dissipative process 

Next, we consider more realistic case with a linear 
damping of the SH wave smaller than that of the FF 
wave, since the FF wave is much closer to the resonance 
where such losses may become quite essential. We as- 
sume again Im(qi^ 2 ) = 0, so that the stationary solution 
of Eqs. ((13|) can be found in this case analytically under 
the conditions of the exact phase matching. First, we 
find the integral of motion of the coupled system in the 
form, 



Vg2 

Vgl 



[516^ + 2^62] =q 2 bj. 



(22) 



Using this invariant, we find the solutions for the ampli- 
tudes of FF and SH waves as follows 



bi(z) = v x 



V g2 



qiq2v g i 



1/2 



csch 



— {z + z ) 

Vgl 



Ql 



2v 1 (z+z )/v 



91 _ 1 



(23) 



In the lossless case, Eqs. (|23|) coincide with Eqs. <|19[1 . 
however, the presence of finite losses changes the type of 
solutions from power to exponential ones. 



A SLAB OF METAMATERIAL 



where 



(Vg2+/3)(v gl -0) 



2 T\ iq2 



1/2 



A 2 = 



0T Ql ' 



(21) 



is the pulse velocity, and T is a characteristic time scale. 
Solutions ((20(1 are valid provided < —v g2 , which means 
that in the moving reference frame both waves should 
be of the same type: either forward or backward. The 
analytical solution ((20|) and l(21|) is characterized by two 
arbitrary constants, and T, so that the corresponding 
family of solutions is two-parametric. 



Next, we study the SHG process for a layer of the 
thickness L (see Fig. 0} and employ a direct numerical 
approach to solve Eqs. Q. First, we rewrite Eqs. Q in 
the dimensionless form 



d 2 H x 



dz 2 
H 2 



+ [e(4(«) - k$] H x = -H 2 H$, 



dz 2 



+ 4 [e(2u)fi(2w) - k 2 x ] H 2 = -QH\, 



(24) 



where the magnetic field is normalized by the value 
u> <j\/ c 2 \a\\ 2 , z is normalized by the value c/oj, and 
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FIG. 6: Transmission coefficient of the FF wave (dashed) 
and reflection coefficient of the SH wave (solid) vs. the slab 
thickness L, for a fixed amplitude of the transmitted wave, 
Hi(L) = 10~ 2 . Right plot shows a blow up of the region 
depicted by a dashed box. 

Q = cri<72/|ci| 2 - We assume that a slab of the left- 
handed material is illuminated by the FF wave with the 
amplitude Hffl , and the SH wave is generated inside the 
slab, so that the reflected and transmitted waves of both 
the frequencies u and 2w appear (see Fig. |3J|, with the 

amplitudes h[ , Hcf , H[ r \ H% , respectively. Solving 
the coupled-mode equations numerically, we present our 
results for the reflection and transmission coefficients de- 
fined as i^, 2a , = HQ/HP and T^, 2uJ = H[%/H^. 

In the calculations presented here we assume that the 
left-handed material is lossless, and we consider the nor- 
mal incidence. We take the following parameters of the 
composite: uq — 2tt x 5 GHz, ui p — 2ir x 7 s -1 , F = 0.3, 
a = 3mm, d = 6mm, U c Rd = 10 5 CGS. For such param- 
eters the exact phase matching takes place at f pm « 5.37 
GHz. Dependences of the coefficients i?2 W , T u and 
vs. the amplitude of the incident FF wave are shown in 
Fig- El O ne can see that the efficiency of the transforma- 
tion of the incident FF wave into the reflected SH wave 
can be rather high. Larger intensities of the incident 



field result in the multistable behavior of the reflection 
and transmission coefficients. 

Figure [S] shows the dependence of the transmission co- 
efficient of the FF wave and the reflection coefficient of 
the SH wave vs. the slab thickness, for a fixed ampli- 
tude of the transmitted wave. We observe multistable 
behavior of the coefficients for thicker slabs. 



VI. CONCLUSIONS 

We have presented a comprehensive study of the spe- 
cific features of the second-harmonic generation in left- 
handed metamaterials with a quadratic nonlinear re- 
sponse. First, we have demonstrated a possibility of the 
exact phase matching between backward propagating FF 
and forward propagating SH waves. Then, we have de- 
veloped an analytical approach based on a novel type of 
the coupled-mode equations in order to characterize the 
process of the harmonic generation from a semi-infinite 
left-handed metamaterial. In particular, we have demon- 
strated that a surface of the left-handed metamaterial 
can operate as an effective mirror that reflects all gen- 
erated SH waves. Finally, we have performed numerical 
simulations of the second-harmonic generation process 
for a finite slab of the metamaterial, and revealed quite 
unusual multi-valued multistable transmissions and re- 
flections. 
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